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ABSTRACT 

Aims. The present work provides the results of the first six years of operation of the systematic night-sky monitoring at ESO-Paranal 
(Chile). 

Methods. The UBVRI night-sky brightness was estimated on about 10,000 VLT-FORSl archival images, obtained on more than 650 
separate nights, distributed over 6 years and covering the descent from maximum to minimum of sunspot cycle n. 23. Additionally, a 
set of about 1,000 low resolution, optical night-sky spectra have been extracted and analyzed. 

Results. The unprecedented database discussed in this paper has led to the detection of a clear seasonal variation of the broad 
band night sky brightness in the VRI passbands, similar to the well known semi-annual oscillation of the Na I D doublet. The 
spectroscopic data demonstrate that this seasonality is common to all spectral features, with the remarkable exception of the OH 
rotational-vibrational bands. A clear dependency on the solar activity is detected in all passbands and it is particularly pronounced in 
the U band, where the sky brightness decreased by ~0.6 mag arcsec"^ from maximum to minimum of solar cycle n. 23. No correlation 
is found between solar activity and the intensity of the Na I D doublet and the OH bands. A strong correlation between the intensity 
of N I 5200A and [OI]6300,6364A is reported here for the first time. The paper addresses also the determination of the correlation 
timescales with solar activity and the possible connection with the flux of charged particles emitted by the Sun. 
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1. Introduction 

Soon after the beginning of VLT science operations in Paranal, 
ESO started an automatic UBVRI sky brightness survey, with 
the aim of both characterizing the site and studying the long 
term trend, in order to detect any possible effects of human ac- 
tivity. The results obtained during the first 18 months of oper- 
ations (Apiil 2000 - September 2001) have been presented and 
discussed in Patat (I2003al hereafter Paper I). This programme, 
which makes use of all scientific images obtained with FORSl, 
is building one of the most extensive, accurate and homoge- 
neous optical sky brightness data sets ever studied. As shown 
in Paper I, these data allow a very detailed analysis, including 
the study of correlations with other parameters, and the investi- 
gation of short, medium and long term variations. The interested 
reader can find exhaustive reviews on this subject in Roach & 
Gordon ( I1973I I and in Leinert et al. ( 119981 ). while references to 
other published sky brightness surveys are given in Paper I. 

Since the publication of the first results, obtained on 174 dif- 
ferent nights close to maximum of solar cycle n. 23, the data 
base has been steadily growing and progressively extending to- 
wards the solar minimum. In this paper I present a global anal- 
ysis run on the whole data set, which includes broadband obser- 
vations taken on 668 separate nights between 20 Apr. 2001 and 
20 Jan. 2007. Additionally, I present and discuss here a set of 
more than 1000 low-resolution long-slit night sky spectra taken 
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with FORSl between May 1999 and Feb. 2005. They were used 
to measure the fluxes of single lines or integrated OH bands. 

The paper is organized as follows. In Sec.|2]l describe the 
observations and the basic data reduction steps for photometry 
and long slit spectroscopy. In Sec.|3]l review the general results 
obtained during dark time, while the correlation with solar ac- 
tivity and the seasonal variations are discussed in Sects. |4]andl5] 
respectively. The spectroscopic analysis is presented in Sects. ^ 
|7]and|8] Finally, in Sec.|9]l discuss the main results and summa- 
rize the conclusions in Sec.fTOl 

2. Observations and Data Reduction 

The data used in this work were obtained with the 
FOcal Reducer/low dispersion Spectrograph (hereafter FORSl), 
mounted at the Cassegrain focus of ESO-Antu/Melipal/Kueyen 
8.2m telescopes (Szeifert 2002). The instrument is equipped 
with a 2048x2048 pixel (px) TK2048EB4-1 backside thinned 
CCD and has two remotely exchangeable collimators, which 
give a projected scale of 0'.'2 and O'.'l per pixel (24/vm x 24/im). 
According to the collimator used, the sky area covered by the 
detector is 6.' 8x6.' 8 and 3.'4x3.'4, respectively. 

2.1. Photometry 

The photometric data set includes 10,432 images obtained in the 
UBVRI passbands with both collimators. The reduction proce- 
dure is described in Paper I, to which I refer the reader for a more 
detailed description, while here I only recap the basic steps. All 
frames are automatically processed by the FORS pipeline, which 
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Table 1. Main properties of the FORS 1 spectral data set. Only 
spectra taken with sun elevation below -18° were included. 



Grism 


Fiher 


Wav. Range 


Resolution 


Dispersion 


N 






(A) 


(A FWHM) 


(Apx-') 




300V 


GG435 


4300-^8900 


12 


2.6 


676 


300V 




3615-8900(*) 


12 


2.6 


163 


600B 


OG590 


3650-6050 


5.3 


1.2 


143 


600R 


GG435 


5390-7530 


4.5 


1.0 


207 



(*) Second order overlapping above 6000 A. 



applies bias and flat-field correction, the latter performed using 
twilight sky flats. Once the instrument signatures are removed 
from the images, the sky background is estimated using the ro- 
bust algorithm described in Patat ( |2003bl) . The photometric cal- 
ibration into the Johnson-Cousins system is then achieved us- 
ing zeropoints and colour terms derived from the observation of 
standard star fields (Landolt |1992| |. regularly obtained as part of 
the FORS calibration plan. Finally, the observed values are cor- 
rected to zenith using the standard procedure (see for example 
Garstang 1 1989 J and logged together with a number of relevant 
parameters. 



2.2. Spectroscopy 

The spectroscopic data set includes a sub-sample of all long-slit 
science data present in the ESO archive whose proprietary pe- 
riod had expired by the time this paper has been written. For the 
sake of simplicity, I have selected only the data obtained with 
the standard resolution collimator and the single-port high-gain 
read-out mode, since this combination is the most used for long 
slit spectroscopy with FORSl. In order to accumulate a wide 
data sample, I have retrieved from the VLT archive all public 
spectra taken with the 300V grism coupled with the order-sorting 
filter GG435 and a long slit I'.'O wide, which is the most fre- 
quently used (676 frames). The wavelength range was extended 
down to about 3600A in the blue by retrieving also all spectra 
taken with the same setup but with no order-sorting filter, for 
a total of 163 frames. To increase the sample, I have retrieved 
also all spectra obtained with the grisms 600B (143 frames) and 
600R (207 frames) coupled with the order-sorting filters OG590 
and GG435, respectively. Also in these two cases the slit was 1" 
. wide. The main characteristics of each setup are summarized 
in Tab. [T] Exposure times range from a few minutes to one hour 

All images were processed within the longslit package of 
IRAl|^ Due to the large amount of data and the purpose of this 
work, the bias subtraction was performed using only a pre-scan 
correction, while flat-fielding effects were neglected. 

Wavelength calibration was achieved using a 2D solution de- 
rived from a set of reference arc exposures. Given the procedure 
adopted for the spectrum extraction, this step is mandatory, since 
in FORS 1 the line curvature can reach a peak-to-peak value of 
about 10 px. If not accounted for, this instrumental feature would 
produce an apparently significant but artificial line broadening 
when collapsing the 2D spectra perpendicularly to the disper- 
sion direction (see next section). 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 
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Fig. 1. Example of FORSl flux calibrated night sky spectrum 
obtained with the procedure outlined in the text, for the 600R 
grism. Main lines and OH bands identifications are given. The 
insert shows an example spectrum obtained with the 300V 
grism. 



2.2.1. Nightglow spectrum extraction and calibration 

After applying the appropriate 2D wavelength solution to all 
frames, the night sky spectrum is extracted. For this purpose I 
have used a robust algorithm to estimate the mode intensity in 
each column perpendicular to the dispersion direction. This im- 
plicitly assumes that most of the pixels are not "contaminated" 
by the contribution of astrophysical objects, which is reasonable 
in the majority of the cases, as verified by direct inspection of the 
whole two-dimensional data sample. This is both a consequence 
of the relatively large slit length featured by FORS 1 (6.' 8 on 2048 
px) and the typical targets observed with this instrument, which 
are very often faint and star-like sources. After visual inspection, 
only a few spectra were removed from the final data set. 

To allow for completely unsupervised line and continuum 
flux measurements, the accuracy of wavelength calibration is a 
mandatory requirement. Possible causes of rigid shifts in the dis- 
persion solution can be identified as instrument interventions, 
turning into movements of the long slit on the focal plane, and 
flexures at large zenith distances. To correct for these problems, 
I have produced a reference night sky spectrum for each of the 
two resolutions I have used, with a typical accuracy (estimated 
on isolated lines) better than 1 A. Then, by means of cross- 
correlation, the zero point of the wavelength scale of each spec- 
trum is automatically corrected at the end of the extraction pro- 
cedure. This ensures that, at this stage, all spectra have maximum 
wavelength errors that do not exceed lA. 

For the absolute flux calibration I have used a set of spec- 
trophotometric standard stars to derive a reference sensitivity 
function s{A), which I have applied to all spectra. Even though 
this does not take into account the changes in sensitivity which 
are mainly due to the aging of reflective surfaces (Patat i2003al l. 
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Fig. 2. Definition of the blue bands Bl, B2, B3 and B4. 

at the wavelengths of interest they are of the order of a few per- 
cent, and therefore can be safely neglected in this context. 

The flux calibration of the extracted spectrum f{A) to physi- 
cal units is finally computed as: 

F{A) - erg s cm A arcsec 

s{A) t w p 

where t is the exposure time (in seconds), w is the slit width 
(in arcsec) and p is the projected pixel scale (arcsec px"'). 

An example of flux calibrated spectra obtained with this pro- 
cedure is presented in Fig. [T] The resulting signal-to-noise ra- 
tio on the pseudo-continuum changes according to the exposure 
time of the original frame, but it is always larger than 100. 

2.2.2. Line/bands flux measurements 

The line flux integration is carried out in a fully automated way 
within a given wavelength window after subtracting the esti- 
mated pseudo-continuum intensity. The integration boundaries 
have been set according to the spectral resolution. For single 
lines, like the [OI]5577, the optimal window semi-amplitude 
has been set to 1.7xFWHM which corresponds to a +4cr from 
the line center. The pseudo-continuum level is estimated in the 
adjacent emission line-free regions. For more complex features 
(Na ID doublet, OH bands), the integration boundaries and the 
continuum region/s have been set ad hoc, and are reported in 
Table |2] 

Following Barbier (1956'), I have also introduced four bands, 
indicated as Bl, B2, B3 and B4, in the blue spectral domain (see 
Table[2|i. The gap between bands B 1 and B2 has been set in order 
to avoid the strong Ca II H&K absorptions (Fig. |2]i. Finally, a 
broad band ranging from 5500A to 7530A has been introduced, 
with the main aim of giving an overall flux estimate. 

As in the case of the broad-band measurements, the line 
fluxes need to be corrected for airmass. Between the two cases, 
however, there is a difference: in fact, while the integrated flux 



Table 2. Integration boundaries and continuum regions used for 
line flux measurements. 



Line/Band 


Line Range 


Cont. Range 




Setup 


NI 5200 


±4o- 


5040-5120 


(CI) 


300V,600B 


OI] 5577 


±4cr 


5480-5520 


(C2) 


300V,600B/R 


NalD 


Dj - 4cr-D| -1- 4iT 


5800-5850 


(C3) 


300V,600B/R 


OI] 6300 


±4(T 


6400-6450 


(C4) 


300V,600R 


01] 6364 


±4o- 


6400-6450 




300V,600R 


OH(6-l) 


6435-6680 


c r\ a f\f\ 

6750-6800 


(C5) 


300V,600R 


OH(7-2) 


6810-7060 


6750-6800 




300V,600R 


OH(8-3) 


7200-7450 


6750-6800 




300V,600R 


OH(6-2) 


8250-8570 


8160-8230 


(C6) 


300V+GG 


02(0-1) 


8605-8695 


8160-8230 




300V+GG 


Bl 


3630-3890 






300V,600B 


B2 


3990-4330 






300V,600B 


B3 


4350-4850 






300V,600B 


B4 


4870-5170 






300V,600B 


Broad Band 


5500-7530 






300V,600R 



within a broad-band filter is the result of extra-terrestrial sources 
(zodiacal light, unresolved stars and galaxies) and emission 
within the atmosphere, in the case of nightglow emission lines 
all the radiation is of terrestrial origin. Practically this coincides 
with setting / = 1 in Eq. C.3 of Paper I, which I have used to 
correct the observed values. 

In general, the signal in the measured features is so high that 
the uncertainty in the line fluxes is by far dominated by the con- 
tamination by unresolved OH lines and, to a smaller extent, to 
the uncertainty on the pseudo-continuum level. 

3. Dark time night sky brightness 

Since the data set includes observations obtained under a wide 
variety of conditions, in order to estimate the zenith sky bright- 
ness during dark time it is necessary to apply some filtering. To 
this aim I have adopted the same criteria described in Paper 
I: photometric conditions, airmass X <1.4, galactic latitude 
\b\ >10°, helio-eclitpic longitude \A - Aq\ >90°, time distance 
from the closest twilight Af,n,, > 1 hour and no moon (fractional 
lunar illumination FLI=0 or moon elevation Iim < -18°). The 
results of this selection, which reduced the number of suitable 
data points to 3736, are summarized in Table [3] As one can see, 
the average values are all within 0. 1 mag from those reported in 
Paper I (see Table 4). In all filters there is a systematic shift to- 
wards darker values, with the only exception of the / band. Since 
the values reported in Paper I were obtained during the sunspot 
maximum and, given the correlation between solar activity and 
night sky brightness shown by Walker ( 1988), Pilachowski et al. 
(1^89); Kiisciunas ( 1990), Leinert et al. ( 1995 ) and Mattila et al. 
( 1996, Krisciunas ( 1997), Krisciunas et al. (2007 ), this behavior 
was indeed expected. Given the time distribution of data points 
(see Fig.|3]l, the average values reported in Table[3]are biased to- 
wards the sunspot maximum phase. For solar activity corrected 
data, see next section. 

Single measurements for the 5 passbands are presented in 
Fig.|4] which shows also surface brightness distributions for the 
dark time (solid line) and global (dashed line) samples. 

4. Sky brightness vs. solar activity 

As first pointed out by Rayleigh d 19281 ) and confirmed later 
on by several other authors (see for instance Rosenberg & 
Zimmerman 119671 Walker 1 19881 Krisciunas 119901: Leinert et 
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Table 3. Zenith corrected average sky brightness during dark 
time at Paranal. Values are expressed in mag arcsec"^. Columns 
3 to 8 show the RMS deviation, minimum and maximum bright- 
ness, number of dark-time data points, expected average contri- 
bution from the zodiacal light, and total number of data points, 
respectively. 



Total 3736 10432 
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Fig. 3. Image distribution along the time interval covered by the 
present work for the global data set (thin line) and for the dark- 
time set (shaded thick line). The vertical dotted line marks the 
extent of the sample presented in Paper I, while the upper scale 
marks the 1st of January of each year. 

al. [19951 Mattila et al.TW; Krisciunas T997', Krisciunas et al. 
120071) many of the emission features in the night sky spectrum 
show a clear dependency on the sunspot cycle. In particular, B 
and V present a peak-to-peak variation of ~0.5 mag arcsec"^ 
during a full solar cycle. Less clear is the behavior at longer 
wavelengths, which are dominated by the OH emissions, whose 
intensity is uncorrelated with solar activity (see Sec.[6]l. 

The data presented here cover the descent from the maxi- 
mum of sunspot cycle n. 23 to the minimum phase, as shown 
in Fig. |5] that displays the monthly averaged Penticton-Ottawa 
solar flux at 2800 MHz (Covington 1969|^ During this interval 
the Solar Flux Density (hereafter S FD) spans from 0.8 to 2.4 
MJy, a range which is very close to that of a full cycle (the solar 
minimum is expected for the end of 2007). Following what has 
been done by other authors (see for instance Leinert et al. 1995), 
I have studied the correlation between the sky brightness nightly 

^ The data are available in digital form at the following web site: 
Ihttp : / /www . drao . nr c. ca/icarus/www/ archive . html 
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Fig. 4. Zenith corrected sky brightness measured at Paranal dur- 
ing dark time (thick dots) from April 1st, 2000 to April 8, 2006. 
The selection criteria are: \b\ >10°, \A - Aq\ >90°, Af„„- >1 hour, 
FLI=0 or hm < -18°. Thin dots indicate all observations (cor- 
rected to zenith). The horizontal dotted lines are positioned at the 
average values of the selected points, while the vertical dashed 
line marks the end of the time interval discussed in Paper I. 
The histograms trace the distribution of selected measurements 
(solid line) and all measurements (dotted line), while the vertical 
dashed lines are placed at the average sky brightness during dark 
time. 



averages and the SFD monthly averages, computed during the 
30 preceding days. The sky brightness measurements have been 
corrected for the zodiacal light contribution computed for each 
data point as in Paper I (Sec. 4) and using the data by Levasseur- 
Regourd & Dumont JTMOl l. 

All passbands show very good linear correlations, an exam- 
ple of which is presented in Fig.|6]for the V filter In order to give 
a quantitative representation of the effect, I have fitted a relation 
of the type m - nio + y SFD to the data. The results are shown in 
Table |4lfor all filters. Besides reporting the zeropoint (mo, mag 
arcsec ), the slope (y, mag arcsec"^ MJy"') and their associated 
statistical errors, the Table includes also the estimated full solar 
cycle variation (Am - |(2.4 - 0.8)7|), the value attained at solar 
minimum, evaluated for SFD-Q.2i MJy (m„„„), the value corre- 
sponding to the average SFD level < SFD >=1.6 MJy (niave), 
the RMS deviation from the best fit relation (cr), the linear cor- 
relation factor (r) and the number of nights used (A^) for each 
filter 

As one can see, the values of Am are smaller than those re- 
ported by other authors: with the only exception of U, which 
reaches about 0.6 mag arcsec"^, all the others show values that 
are smaller than 0.4 mag arcsec"^. Walker ( I1988I I quoted max- 
imum ranges of AV ^1.0 and AB ^0.8 mag arcsec"^ for solar 
cycle n. 21, while Krisciunas (119971) reports Ay=0.6 for solar cy- 
cle n. 22, and similar values are reported by Leinert et al. ( 1993) 
and Mattila et al. (1996). On the other hand, Liu et al. (2001} 
quote an increase of the V sky brightness of ~0.2 mag arcsec"^ 



Filter Sky Br. cr Min Max Nj, 



U 
B 
V 
R 
I 



Amz 



22.35 0.19 

22.67 0.16 

21.71 0.24 

20.93 0.24 

19.65 0.28 



21.89 
22.19 
21.02 
20.42 
18.85 



22.78 
23.02 
22.30 
21.56 
20.56 



129 
493 
692 
1285 
1137 



0.20 
0.28 
0.20 
0.16 
0.07 



Nra 



264 
1400 
1836 
3931 
3001 



22 
21 



22 
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Fig. 5. Penticton-Ottawa Solar flux at 2800 MHz (monthly aver- 
age). The time range covered by the data presented in this paper 
is indicated by the thick line, while the horizontal segment in- 
dicates the time covered by the data presented in Paper I. The 
upper insert traces the solar flux during the last six cycles. 
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Fig. 6. Lower panel: nightly average dark time sky brightness in 
the V passband vs. solar flux density. The solid line is a linear 
least squares fit to the data (empty symbols mark the data points 
rejected by a kappa-sigma clipping). Upper panel: Penticton- 
Ottawa solar flux at 2800 MHz during the time interval discussed 
in this paper. The open circles indicate the monthly averaged val- 
ues corresponding to the V nightly averages plotted in the lower 
panel. 



Fig. 7. Upper panel: linear correlation coeflicient (absolute 
value) as a function of time delay r for the V passband data, 
computed for three different values of At (1, 10 and 30 days). 
Lower panel: example of SFD periodicity due to solar rotation. 
The dotted curve is a replica of the original data shifted by 27.3 
days. 



from 1995 to 2001. This value is consistent with the measures 
discussed here, especially taking into account that cycle n. 23 
had a second maximum, which occurred after the observations 
presented by Liu et al. (2003). These facts seem to suggest that 
not all solar cycles have identical effects on the night glow. As a 
matter of fact. Walker ( 1988), while revising the result of previ- 
ous works, had suggested that the relation between intensity of 
the [OI]5577 A line and the solar activity might vary from cy- 
cle to cycle, within a given cycle and possibly with geographical 
location. Unfortunately, the number of sunspot cycles covered 
by the observations is still too small to allow a firm conclusion, 
but the very recent results discussed by Krisciunas et al. ( 20071 1. 
covering two full solar cycles, seem indeed to confirm this sus- 
picion. 

4. 1. Solar activity correlation time scales 

Given the unprecedentedly large amount of data, one can investi- 
gate the solar dependency in a bit more detail, for example trying 
to deduce the typical timescales of night sky brightness fluctu- 
ations due to short-term changes in the solar flux. Ideally, for 
doing this, one would look for strong enhancements in the solar 
flux and try to seek for a corresponding increases in the night 
sky brightness. Nevertheless, due to the sparse time sampling, 
the only viable approach is actually the opposite one, i.e. start- 
ing from the available night sky measurements, one goes back 
to the solar data and studies the correlation with the sun flux 
emitted as a function of time lag. For doing this I have com- 
puted the Pearson's linear correlation coefficient r (Press et al. 
J992) between the nightly average sky brightness m(t) measured 
at any given time t and the solar flux density SFD{t - t), aver- 
aged within a time window At, as a function of both t and At. 
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Table 4. Linear least squares fit parameters for the sky bright- 
ness vs. solar activity relation. Values in parenthesis indicate 
RMS uncertainties. Input data have been corrected for the dif- 
ferential zodiacal light contribution (see text). 



Filter m„, 



Am 



mo 



N 



U 22.86 22.58 0.61 23.15 (0.12) -0.36(0.07) 0.15 0.47 32 

B 23.11 22.98 0.29 23.25 (0.04) -0.17 (0.02) 0.12 0.40 127 

V 21.99 21.86 0.30 22.13 (0.05) -0.17 (0.03) 0.14 0.42 148 

R 21.26 21.09 0.37 21.33 (0.04) -0.22(0.02) 0.15 0.44 202 

I 19.81 19.72 0.20 19.90(0.06) -0.11 (0.04) 0.18 0.28 144 



While the first parameter gives an indication about the time lag 
between a change in the solar flux and the consequent variation 
in the night sky brightness, the second is related to the typical 
timescales of the physical processes which govern the energy 
release in the upper atmosphere. 

An example of this kind of analysis is shown in the upper 
panel of Fig. [7] which illustrates the behavior of the correlation 
coefficient r for the V data, obtained on 170 separate nights for 
three dilTerent values of At (1, 10 and 30 days). The correlation 
peak is quite broad and it is centered at r ~30 days; moreover, 
the maximum correlation tends to increase slightly for larger av- 
eraging windows. This plot also explains why the average solar 
flux computed in the 30 days before the night sky observations 
(empty circles) gives a better correlation than the sun flux mea- 
sured on the preceding day (filled circles), as already pointed out 
by several authors (see for example Leinert et al. 119951 Mattila et 
al. 1 19961 ). Finally, Fig.[7]clearly shows that the correlation drops 
significantly for t >40 days, a fact that is common to all UBVRI 
passbands. 

An interesting feature to be noticed in Fig.[7](seen also in the 
equivalent plots for the other filters), is the presence of spurious 
correlation peaks at a constant separation of about 27 days from 
the main peak. This is due the a periodicity present in the solar 
flux data (see Fig. [7] lower panel), which is related to the solar 
rotation, whose synodic period is ~27.3 days (Howard I1999I I. 
This semi-regular recurrence in the solar data explains, for ex- 
ample, the presence of the two bumps close to t-5 and t-65 in 
Fig.|Z] 

The strongest correlation is shown by the U passband data, 
which presents a rather marked peak r ~0.6 at t ~15 days, while 
r ~0. 15 for t=1 day. For this reason, averaging over the last 
30 days gives a strong increase in the correlation, much stronger 
than in any other passband. A behaviour similar to that displayed 
in V is seen also in B (correlation peak r ~0.45 for t ~25 days) 
and / (r ~0.30 for t ~20 days). Somewhat different is the case 
of R passband, for which the correlation peak (r ~0.45) is at- 
tained at T ~2 days, suggesting that the sun-dependent features 
that contribute to the flux in this filter react rather rapidly to the 
solar flux fluctuations. In general, however, and with the possi- 
ble exception of the U passband, the correlation peaks are rather 
broad, indicating that different processes take place with differ- 
ent timescales. The behavior of the U band, where the night- 
glow emission is dominated by the Herzberg and Chamberlain 
O2 bands (Broadfoot & Kendall [T968|l, indicates that the photo- 
chemical reactions that are responsible for the emission in this 
region are more sensitive to solar activity. 



5. Seasonal Variations 

In the previous work I had attempted to detect night sky bright- 
ness seasonal variations but, due to insufficient number of data 



points, I could not draw any firm conclusion (see Fig. 14 in Paper 
I). Thanks to the much larger sample now available, this analy- 
sis becomes feasible and, as a matter of fact, traces of a periodic 
modulation in the average sky brightness are visible already in 
Fig.|4] They become much clearer when each data point is plot- 
ted against the number of days from the beginning of the cor- 
responding year The result is shown in Fig. [8] where, besides 
reporting the single dark time measurements, I have also plotted 
the monthly averages. The input data have been corrected for 
differential zodiacal light contribution and the solar flux depen- 
dency derived in the previous section has been removed using 
the parameters presented in Table [4] 

This semi-annual oscillation (hereafter SAO) is definitely 
present in V, R and /, while its presence in B is more question- 
able (U data were not included since the sample in this passband 
is too poor for this purpose). The modulation amplitude grows at 
longer wavelengths, shows two maxima around April-May and 
October, and two minima around July-August and December- 
January. In general, the variation is more pronounced in Winter- 
Spring than in Summer-Fall. For example, in / it reaches a peak- 
to-peak value of about 0.5 mag arcsec"^. 

As pointed out by Benn & Ellison ( 19981 1. the variable contri- 
bution of zodiacal light can mimic a seasonal variation. In order 
to exclude a possible contribution by this source to the observed 
behavior, I have analyzed the expected enhancement of bright- 
ness due to the zodiacal light for each data point, using the data 
presented by Levasseur-Regourd & Dumont (1 19801 1 and the pro- 
cedure discussed in Paper I (Sec. 4) for the Paranal site. This 
does not show any significant regular pattern as a function of 
day of the year The conclusion is that the observed SAO is not 
due to the periodic apparent variation of the ecliptic height above 
Paranal's horizon. 

An interesting thing to be noticed, is that the minima and 
maxima of the SAO occur out of phase with respect to the 
Equinoxes and Solstices (see Fig.|8] vertical dotted lines). 

While seasonal variations of emission lines and/or bands 
have been studied by several authors in the past (see for exam- 
ple Chamberlain 1961 Roach & Gordon 1973' and references 
therein), broad band measurements are much more scanty and 
the results not always in agreement. For instance, Schneeberger, 
Worden & Beckers (119791 ) report particularly bright values ob- 
tained in June at the Sacramento Peak Observatory and they find 
them to be marginally correlated with the strong seasonal trend 
evident in the record of daytime sky brightness observations. In 
their survey run at the Lowell Observatory, Lockwood, Floyd & 
Thompson (119901 ) discuss the seasonal variation, concluding that 
neither winter enhancements [...] nor springtime rise [...] is in- 
dicated [...]. Benn & Ellison ( 119981 ) reach the same result from 
the analysis of the data obtained on La Palma, concluding that 
dark-of-moon sky brightness does not vary significantly (<0.1 
mag) with season [...]. Finally, Liu et al. ("2003"), analyzing data 
taken at the Xinglong Station between 1995 and 2001, find that 
the sky is darker in the fall and winter than in the spring and 
summer [...]. While part of the discrepancies can be due to lat- 
itude effects (see for example Chamberlain 119611) . some of the 
negative detections are probably to be ascribed to non sufficient 
time sampling and coverage. In fact, the SAO amplitude is at 
most comparable with the night-to-night fluctuations and hence 
large and well sampled data sets are required. 

It is worth mentioning here that Garstang (119881 ) has pro- 
duced some simplified models to predict seasonal variations in 
the broad-band night sky brightness, based on periodic varia- 
tions of height and molecular density. The values predicted by 
these models are far too small with respect to those presented 
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Fig. 8. Seasonal variation of dark time sky brightness with re- 
spect to the average value. Small symbols are the single mea- 
surements while the large symbols mark the monthly averages, 
computed in bins of 30 days each. The error bars indicate the 
RMS deviations from the average within each bin, while the ver- 
tical dotted lines mark equinoxes and solstices. Data have been 
corrected for differential zodiacal light contribution and solar 
flux dependency. 
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Fig. 9. Lower panel: dark time, zenith corrected [OI]5577 line 
flux as a function of S FD. The solid line traces a linear least- 
squares fit to the data. Upper panel: line fluxes as a function 
of time from the beginning of the year The data have been 
corrected to solar minimum {SFD-Q.% MJy) using the relation 
shown in the lower panel. The large points mark monthly aver- 
ages and the right vertical scale is expressed in Rayleigh. 



here, suggesting that other possible explanations must be inves- 
tigated. 

The observed behavior might indicate that whatever the rea- 
son for the periodic variation is, it is not directly related to the 
amount of sun radiation received by a given patch of the atmo- 
sphere during the day. In fact, one might think that since dur- 
ing the austral summer days are much longer than nights, this 
could result into a brighter nightglow. The data show actually 
the opposite behavior, since during austral summer the night sky 
reaches its lowest average brightness. Moreover, this appears to 
be in phase with what is observed in the northern hemisphere, 
where the sky is darker in winter than is summer (see for exam- 
ple Liu et al. 2003 ). This seems to indicate that the SAO must 
be related to some other, non local mechanism (see the discus- 
sion in Sec. |9]). Remarkably, but after all not surprisingly, polar 
auroral activity shows a similar temporal fluctuation, with max- 
ima in spring and autumn (see for example Meinel, Neighed & 
Chamberlain [19541 ■ 

6. Spectroscopic analysis 

Since there are many distinct components that contribute to the 
global nightglow emission (Leinert et al. 1998 and references 
therein) that cannot be disentangled with broad-band photome- 
try, several attempts have been made in the past to increase the 
spectral resolution using a set of narrow band filters. After the pi- 
oneering work by Barbier ( 1956; see also Chamberlain '1961 for 
a review on this subject), who used 8 intermediate-band photom- 
etry, several other researchers have undertaken similar projects 
(see for example Leinert et al. 1 1995 1 and Mattila et al. 1 19961 for 



two more recent works). Nevertheless, due to the relatively large 
bandwidths (typically broader than 100 A), the study of single 
features, especially if not very intense, has always been hindered 
by the presence of the pseudo-continuum and possible adjacent 
lines and/or bands. 

6.1. Main atomic emission lines 

In this section I will analyze the behavior of the main emission 
lines in the optical domain. For convenience, I will express the 
line fluxes in Rayleigh^ 

6.1.1. The [01] lines 

The [OI]5577 is generally the most prominent feature in the op- 
tical night sky spectrum. It falls right in the center of the V pass- 
band, giving a typical contribution of 20% to the global surface 
brightness in this filter. It has a typical intensity of ~250 R, it 
arises in layers placed at about 90 km (Roach & Gordon I 1973I I 
and it displays a marked dependency on solar activity (Rayleigh 
11928 ). This is clearly shown also by the data presented here 
(Fig. |9] lower panel), which indicate also the presence of pro- 
nounced fluctuations (40 to 750 R peak-to-peak) around the av- 
erage level (~230 R). To quantify the correlation with solar ac- 
tivity and following the procedure that has been applied to the 
broad band data (see Sec. [4|, I have fitted to the data a law of 
the type log F = log Fq + y S FD. The results for this and other 
lines are presented in Table |6j which shows also the number of 

3 1R=10V4:t photons s"' cm-^ sr"' = 3.72x10-'"' /i"'(A) erg s"' 
cm"^ arcsec"^. 
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Table 5. Observed line intensities and relative flux contribution 
to V and R passbands for [OI]5577, [OI]6300 and Na I D. All 
values are corrected to zenith. 



[OI]5577 
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(10"'* erg s"' cm"^ arcsec"^) 
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40 
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<1% 


ave 
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230 
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1% 


max 


50 
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62% 


4% 



Table 6. Linear least squares fit parameters for the line fluxes 
vs. solar activity relation. Values in parenthesis indicate RMS un- 
certainties. Input data have been corrected to zenith. 



Line logFo y a r N 

OI 5577 -15.36 (0.03) 0.27 (0.02) 022 0:49 876 

01 6300 -15.69(0.07) 0.25 (0.05) 0.41 0.23 635 

NalD -15.72(0.03) 0.04(0.02) 0.25 0.07 876 

NI5200 -17.19(0.05) 0.37 (0.03) 0.32 0.45 727 
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data points used (A^) and the linear correlation coefficient (r). 
What is interesting to note is that, once the solar dependency is 
removed from the data, the flux of this line displays a marked 
SAO (see Fig. [9) upper panel), similar to that seen in the broad 
band data (FigTSll. This fact has been already noticed by Buriti 
et al. ( 1200 11 1, who found that this line and other mesospheric fea- 
tures all show a SAO. 

The same behavior is, in fact, shown by the [OI]6300,6364A 
doublet (see Fig. 10 1, which is produced at 250-300 km (Roach 
& Gordon |1973l l and it is known to undergo abrupt intensity 
changes on two active regions about 20° on either side of the 
geomagnetic equator (Barbier [l957l l. hence marginally includ- 
ing the Paranal site. Indeed, the [OI]6300 measured fluxes show 
very strong variations (10 to 950 R) around the average level 
(~150 R), with spikes reaching ~1 kR, so that the line fluxes 
span almost a factor 100 (to be compared with the factor ~25 
measured for the [OI]5577. In this respect, the behavior of the 
[OI]6300,6364A doublet is different from that of the [OI]5577 
line, since in a significant number of cases its flux is very small 
and the line is practically invisible, lost in the OH(9-3) molecular 
band (see Fig.[T]i. In fact, its flux is less than 80 R for more than 
50% of the cases considered in this work, the minimum recorded 
value being ~10 R (Fig. [TT) . The variation range of these two 
features and their relative contributions to the V and R total flux 
are summarized in Table |5] From the values reported in this ta- 
ble, one can see that [OI]5577 can produce a maximum variation 
of ~ 0.5 mag arcsec"^ in the V band, while a similar effect is pro- 
duced by the [OI] doublet in the R passbancQ Given the fact that 
the RMS variation in the V and R passbands is about 0.25 mag 
arcsec"^ (see Table |3]l, this implies that the fluctuations seen in 
these passband are not completely accounted by the changes in 
the atomic O line fluxes. The two lines appear to show a very 
weak correlation: the average ratio F(6300)/F(5577) is 0.64, but 
in a significant number of cases (~25%) this ratio is larger than 
1 (see also Sec. [tJi. In those circumstances, the [OI]6364A line 
becomes the most prominent nightglow atomic feature in the op- 
tical domain. 

As far as the solar activity is concerned, it must be noticed 
that even though a dependence from the SFD is seen, it is less 
clear than in the case of the [OI]5577 (see Fig. 



10 lower panel). 



In fact, even though the slopes are very similar, the linear corre- 



* Since the [OI]6364/[OI]6300 lines ratio is 1/3, the maximum com- 
bined contribution is about 60%. 
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Fig. 10. Same as Fig.^for the [OI]6300 fine. 



lation coefficient for [OI]6300 is 0.23, to be compared with 0.49 
measured for the [OI]5577 line (see also Tab.|6]l. 

It is well known that, during the first hours of the night, the 
[OI]6300 line shows a steady decline in brightness which, at 
tropical sites, is interrupted by abrupt emission bursts (see, for 
instance. Roach & Gordon 1 19731 . A plot of the line fluxes as a 
function of time elapsed after the evening twilight (see Fig. 12 1 



clearly shows that these events occur during the whole night, and 
not only before midnight, as sometimes stated in the literature 
(see for instance Benn & Ellison IT998I ). Actually, one interest- 
ing fact that emerges from this analysis is that, when the nights 
are short (filled symbols), i.e. during austral summer, most of 
the line enhancements tend to take place in the second half of 
the night. On the other hand, when the nights are long (i.e. last- 
ing more than 8.5 hours), flux changes tend to become smaller 
and smaller during the last two hours of the night. No such di- 
chotomy is observed for the [OI]5577 line (see Fig. 
panel). 



12 lower 



Finally, as for the broad band data (see Sec. 4. 1 1, 1 have cal- 
culated the correlation coefficient between line fluxes and SFD 
for different values of the time delay t and time window At. 
An instructive example is presented in Fig. 13 showing the case 



of the [OI]5577 line, which displays the strongest dependency 
on solar activity. The correlation function shows a peak at about 
T=15 days for At=1 day. The maximum correlation increases for 
a At=10 days, peaking at t ^12.5 days. As for the photometric 
data (see Fig. |7]i, spurious correlation peaks due to the 27 days 



F. Patat: The Dancing Sky: 6 years of night sky observations at Cerro Paranal 



9 




Fig. 11. Dark time line flux distributions for [OI]5577 
[OI]6300,6364A, Na I D and N I 5200. The vertical lines trace 
the average (dashed) and the median (dotted-dashed) of the dis- 
tributions. The upper horizontal scale in each panel is expressed 
in Rayleigh. 



Fig. 12. Line fluxes as a function of time elapsed from evening 
twilight for [OI]6300 (upper panel) and [OIJ5577 (lower panel). 
The vertical dashed lines mark the shortest and longest nights 
(7.4 and 10.7 hours respectively), while dotted lines indicate the 
midnight in the two cases. Filled symbols indicate data obtained 
in nights whose length is between 7.4 and 8.4 hours. 



solar rotation are present. A similar analysis for the [OI]6300 
line shows two similar peaks at r ^13 days and 27 days later. 

6.1.2. The Na I D lines 

The Na I D doublet is known to originate in a layer placed at 
about 92 km and to undergo strong seasonal variations around an 
average value of 50 R (Roach & Gordon 1973). For the Paranal 
site, the doublet varies by a factor larger than 15, reaching a max- 
imum value of 160 R in April. The maximum contribution of the 
Na I D doublet to the global night sky brightness reaches ~0. 1 
mag arcsec"^ both in V and R passbands (see Tablejsjl. 

While the Na I D doublet intensity does not show any signif- 
icant correlation with solar activity (the correlation coefficient is 
only 0.07; see Fig. 14 and Table |6|l, it does show a clear SAO, 
with a peak in April and a secondary peak in October/November 
(see Fig. [T4j upper panel), as it is seen in all atomic lines ana- 
lyzed in this work. The SAO displayed by mesospheric sodium 
column density is a well studied phenomenon (Kirchoff 1 19861 1. 
also because of the importance of the sodium layer for the laser 
guide star adaptive optics systems (see for instance Ageorges & 
Hubin BOOOl l. The fact that different species like Na I and [OI] 
all show a SAO is interpreted in the light of the role of O atoms 
in the source photochemical reactions (Slanger 2005 ). 

Since the resolution provided by the grism 600R is sufficient 
to resolve the Di and D2 components (see Fig. 15 1, following 
the work done by Slanger et al. (2005 ), I have estimated the in- 
tensity ratio D2/D1 during dark time, on a total of 147 spectra. 
The region of interest is contaminated by the presence of at two 
features belonging to the OH(8-2) band, which I will indicate 
as B and C (see Fig. 15 for their identification). Using a high 



resolution UVES spectrum of the night sky (Hanuschik 2003), I 
have estimated the intensity ratio between these two features and 



another OH(8-2) line at 5932.9A (Pii(3.5), Abrams et al. [T994] l, 
which I will indicate as E. Since this feature is well measurable 
in the 600R spectra, assuming that these intensity ratios are con- 
stant (B/E=0.95, C/B=0.35), I have estimated the intensity of B 
and C from the measured intensity of component E. Then, after 
subtracting to the original data two Gaussian profiles centered at 
the laboratory wavelengths of B and C, I have fitted the resid- 
ual with a double Gaussian profile, fixing the distance between 
the two components (5.98A) and their FWHM (4.5A). FinaUy, 
I have derived the D2/D1 simply computing the ratio between 
the central intensities of the two fitted Gaussians. The pseudo- 
continuum has been removed using a first order polynomial. 

The results, shown in Fig. [16] (upper panel), are perfectly 
in line with the findings published by Slanger et al. (i2005) : the 
intensity ratio, which should be 2.0 if the two transitions are 
produced according to their spin-orbit statistical weights, varies 
from 1.2 to 1.8, with most of the data lying between 1.5 and 1.7. 
The average value for the FORSl sample is 1.64, with an RMS 
deviation of 0.08. Using a larger sample including more than 
300 high resolution spectra, Slanger et al. (12005 l l have found 
evidences for a SAO for the D2/D1 intensity ratio. This is not 
detected in the low-resolution data set presented here, but this 
is probably due to the fact that the present sample includes less 
than half as many data points. 

6.1.3. The N I 5200 line 

The N I feature at ~5200A is actually a blend of several tran- 
sitions. It is supposed to originate at about 260 km altitude and 
has a typical intensity of 1 R (Roach & Gordon 1973 ). As I have 
shown in Paper I, this line shows abrupt changes, possibly fol- 
lowing the behavior of the [OI]6300 line. The data discussed 
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Fig. 13. Linear correlation coefficient (absolute value) as a func- 
tion of time delay t for the [01]5577 line, computed for four 
different values of At (1, 10, 30 days). 



JFMAMJ JASOND 



-15.0 



K -16.0 - 




E 

o 

-15.0 " 



1 00 200 300 

Days from January 1st 



-16.5 



— I 1 r 




_i I I I I L_ 



2 2 
Solar Flux Density (MJy) 



Fig. 14. Same as Fig. [9] for the Na I D doublet. In this case, no 
solar flux correction has been applied. 



in this paper show that the flux of this line ranges from prac- 
tically zero (the line is not detected) to about 30 R; moreover, 
the flux distribution is rather similar to that of [Ol]6300 (see 



Fig. Ill, strengthening the impression that these two lines are 



related. This line displays also a strong dependency on solar ac- 
tivity, with a correlation factor similar to that of the [01]5577 
line (see Table|6]l. Finally, as for all other atomic lines discussed 
here, it shows a SAO (see Fig.fTT]). 



Fig. 15. Night sky spectrum taken with UVES (Hanuschik l2003l l 
in the region of the Na 1 D lines (resolution 0.15A FWHM. 
The labels A, B and C indicate the OH Meinel 8-2 band lines 
222(0.5) 5887 lA, 6ii(1.5) 5888.2A and eii(2.5) 5894.5A 
(Abrams et al. 119941) . The insert shows the corresponding spec- 
trum obtained with FORS 1 and the 600R grism (resolution 4.5 A 
FWHM). The same spectrum is shown also in the main panel 
(thin line). 
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Fig. 16. Upper panel: distribution of the Na I D2/D1 intensity 
ratio during dark time. Lower panel: Na 1 D2/D i intensity ratio as 
a function of time from the beginning of the year The horizontal 
dashed line marks the average value. 
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Fig. 17. Same as Fig.^for the N I 5200A line. 



Fig. 18. Same as Fig.|9]for the OH(8-3) molecular band. No solar 
flux correction has been applied. 



The maximum contribution of this line to the global flux in 
the V band during dark time is about 2%. 

6.2. Main molecular emission bands 
6.2.1. The OH bands 

The OH bands analyzed here (see Table|2| show a very tight mu- 
tual correlation, in the sense that they appear to vary in unison. 
Moreover, they do not show any correlation with solar activity 
(r <0. 1 for all bands), while they show a SAO, even though not 
as pronounced as in the case of the other features discussed so far 
(see Fig. 18 for an example). The flux distribution appears to be 
much more symmetric around the average value than in the case 
of atomic lines. All bands shows the same range of variation, 
which is close to a factor 2 around the mean value. Given the in- 
tensity of these features (especially OH(8-3) and OH(6-2)), their 
variability is certainly the dominating source of sky brightness 
fluctuations in the I passband. 



6.2.2. The 02(0-1) band 

This band shows a clear correlation with solar activity (r=0.41) 
and the same SAO observed for all other features (see Fig. 19 1. 
Its integrated flux varies by about a factor 15, between 80 R and 
1.1 kR. 
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Fig. 19. Same as Fig.|9]for the 02(0-1) molecular band. 



6.3. The blue bands and pseudo-continuum 

None of the blue bands, which are supposed to trace the be- 
havior of the pseudo-continuum generated by the O2 and NO2 
molecular bands (Roach & Gordon fl 9731 1. shows a clear depen- 
dency on solar activity (r <0.3). However, it must be noticed 
that, while the emission features are purely generated within 
the atmosphere, the continua are significantly influenced by the 
extra-terrestrial background (Roach & Gordon 1973), which is 
difficult to remove. In fact, for the sake of simplicity, no differ- 



ential zodiacal light contribution has been applied to the con- 
tinuum measurements. This is certainly affecting the blue bands 
and the continuum regions, especially those close to the zodiacal 
light spectrum peak. This and the improper removal of airmass 
effects is most likely the cause for the lower correlations shown 
by these features and the solar activity, with respect to what is 
detected, for example, for the U and B passbands. 

The integrated flux varies of about a factor 2 around the av- 
erage value, with a fairly symmetric distribution. As already in- 
dicated by the B band data (see Sec. |5]l, there is no significant 
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Fig. 20. Seasonal variation of blue band B4 (upper panel) and 
continuum range C2 (lower panel). 



trace of a SAO for none of the blue bands (see Fig. 20 lower 
panel). Interestingly, while also the continuum range CI (5040- 
5120A, see Table does not display any evidence for seasonal 
fluctuations, the redder ranges C2 to C6 show a possible broad 
annual oscillation, with a peak in June (Fig. 20 lower panel). 

While the region between 3100 and 3900A (mostly not cov- 
ered by the data set discussed in this paper), shows a number of 
emission features attributable to O2 Herzberg and Chamberlain 
bands, the spectral interval 3900-4900 is almost a pure pseudo- 
continuum (see for instance Broadfoot & Kendall 119681 and 
Fig.|2]here). Nevertheless, clear variations are seen in the FORS 1 
database, as it is illustrated in Fig. |2T| where 1 have compared 
to spectra obtained with the same instrumental setup and simi- 
lar signal-to-noise ratio on two different nights. Whilst the two 
spectra show the same overall emission features, these are much 
more pronounced in the data obtained in November 1999 . 

For an explorative analysis, 1 have measured the fluxes of 
some of the most prominent features, which 1 have marked in 
Fig. 



21 



(lower panel) and I will indicate as bfl (3720-3767 A), 
bf2 (3807-3860A), bf3 (4150-4200A), bf4 (4387-4460A) and 
bf5 (4809-4859 A). While bfl, bf2 and bf3 can be tenta tively 
identified as N2 and N:J^ molecular bands (see Chamberlain 11961] 
and also the next section here), the identification of bf4 and bf4 
is more difficult. For simplicity, 1 have used a common value 
for the blue pseudo-continuum (CB), which was estimated in 
the spectral range 4730-4760A, that does not show any trace of 
emission features (see for instance Hanuschik l2003l l. The results 
are summarized in Table |7] that reports the maximum and me- 
dian values derived from the 114 suitable spectra used in the 
analysis, together with the percentage of non-detections. 

The largest range of variation is shown by bfl, which can 
reach a maximum flux of about 1.4x10"'^ erg s"' cm"^ arcsec"^, 
while bf2 to bf5 reach peak fluxes between 6 and 8x10"'^ 
erg s"' cm"^ arcsec"^. Despite their similar peak values, these 
features display different behaviors. For instance, while bf5 is 
practically always detected, bf2 is absent in more than 50% of 



Fig. 21. The spectral region 3750-5400A on two different 
epochs. Both spectra were obtained in dark time with the 600B 
grism. The two narrow absorption visible between 3900 and 
4000A are the Call H&K solar features. 



Table 7. Dark time maximum and median fluxes for the blue 
emission features indicated in Fig. 



21 The last column gives 



the frequency of non-detection (flux below 10 
arcsec"'). 



erg 



bf 


wavelength range 


max 


med 


n.d. 


# 


(A) 


(10"" erg s" 


cm"^ arcsec"^ ) 




1 


3720-3767 


14.4 


5.4 


26% 


2 


3807-3860 


8.4 


0.0 


54% 


3 


4150-4200 


6.4 


1.3 


38% 


4 


4387-4460 


7.2 


2.8 


17% 


5 


4809-4859 


8.0 


3.9 


1% 



the cases. Since the median value of the pseudo-continuum is 
7.2x10"'^ erg s"' cm"^ arcsec"', each of these features con- 
tributes less than 2% to the integrated flux between 3700 and 
5000A. Therefore, their fluctuations can account only partially 
for the variations observed in the B passband, that must be re- 
lated to the changes in the pseudo-continuum. 



6.3.1 . The strange event of November 8, 2004 

Normally, there are no conspicuous isolated emission features 
in the airglow bluewards of 5200A, which is dominated by 
the pseudo-continuum (Roach & Gordon 1973 See also Fig. [2] 
here). However, N. Castro and M. Garcia, while analyzing a set 
of low resolution spectroscopic data obtained with FORS2 on 
November 8, 2004, have noticed the presence of unexplained 
emission features between 3600 and 4400A. A more detailed 
analysis has shown that these emissions where present on all the 
FORS2 data of that night, i.e. a set of MXU exposures 2700 
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Fig. 22. The weird night sky spectrum of November 8, 2004. 
The first negative band identifications are from Chamberlain 
( 119611 Table 5.4). For presentation the Nov 11 spectrum has 
been scaled by a factor 1.1, in order to match the pseudo- 
continuum on Nov 8. 



seconds each, obtained between 00:15 and 07:08 Ulj^using the 
600B grism. As expected, no trace of these features was visi- 
ble in a similar data set obtained three days later An example 
is shown in Fig. 22 where, for comparison, an analogous spec- 



trum obtained on November 1 1 is also plotted. Clearly, the two 
spectra diff'er mainly for the presence of two prominent emission 
bands, peaking at 3194A and 4278A, which are identified asNJ 
first negative bands IN(O-O) and IN(O-l) (Chamberlain [T96il 
Table 5.4). These features, which are normally very weak or 
even absent in the nightglow (Broadfoot & Kendall [19681 1 are on 
the contrary typical of the aurora spectrum (Chamberlain 1 1 96 il l . 
Besides being an extremely strange phenomenon at the latitudes 
of Paranal, an aurora would certainly be accompanied by other 
spectral markers, like for instance a large increase in the emis- 
sion of the [OI]5577 line, which can reach in fact an intensity 
of 100 kR during a IBC III aurora (Chamberlain irMTT) . The flux 
carried by this line in the same spectrum presented in Fig. |22] 
is ~183 R, that is slightly below the average level measured for 
Paranal (230R, see Sec. 



6.1.1 



. This definitely rules out an excep- 
tional auroral event as the responsible for the unusual spectrum 
observed on November 8, 2004. 

A plausible explanation, proposed by Castro & Garcia (pri- 
vate communication), is the contamination by the reflection from 
clouds of a number of lightning strokes. In effects, in the spec- 
tral range covered by the F0RS2 data (3600-6 100 A), the most 
prominent features in a lightning spectrum are the Nj first nega- 
tive bands IN(O-O) and IN(O-l) (see for instance Wanace |1964| l. 
Additionally, on the night of Nov 8 2004 thick and thin cirrus 
were reported in the ESO-Paranal night logs, substantiating the 
hypothesis of scattered light from a rather far thunderstorm. 



^ On November 8, 2004 the evening twilight ended on 00:23 UT and 
the morning twihght started on 08:27 UT. 



This kind of events must be indeed very rare, since no other 
example could be found in the FORSl spectral data base pre- 
sented in this paper. 



7. Correlations between spectral features 

For a first exploratory analysis I have computed the linear corre- 
lation coefficients in the log F-log F plane between all measured 
features.. The results are presented in Table [8] and they basically 
confirm the correlations found by Barbier (|1956l l, even though 
new interesting facts do appear 

The pioneering optical, eight-color photometric studies by 
Barbier (1956) have shown the existence of the so-called co- 
variance groups: the green-line group ([OI]5577, O2 Herzberg 
bands, the blue bands, the green continuum and the ©2(0-1) 
band), the sodium group (Na I D doublet and the OH bands) 
and the red-line group, which includes only the [OI]6300,6364 
doublet). So far, the latter appeared to be completely indepen- 
dent from any other component of the airglow (see Chamberlain 
119611) . Nevertheless, as it is shown in Fig. |23] [OIJ6300 shows 
a very tight correlation with the N I feature at 5200A. The lin- 
ear correlation factor in the log-log plane is r=0.95 and this ap- 
pears to be one of the strongest correlation between airglow fea- 
tures found in the data set presented here, surpassed only by that 
shown by the OH bands (see Table [8]l. To my knowledge, this 
is the first time this finding is reported; most likely, it escaped 
the attention of previous investigations simply because the N I 
feature is rather weak (<30 R) and hence practically impossible 
to measure with intermediate passband filters. 

Even though a correlation between N 1 5200 and [OI]5577 is 
found (see Fig. |23] upper panel), this is less marked (r=0.56) 
and the spread around the best fit relation is much larger 
(cr=0.29 vs. cr=0.1 1). Finally, the correlation between [OI]5577 
and [OI]6300 is indeed weak (r=0.29, cr=0.39); neverthe- 
less, the data presented here seem to indicate that, on av- 
erage, the maximum value attained by the red line is re- 
lated to the flux of the green-line through the simple relation 
F([OI]6300)<2xF([OI]5577). 

The strongest correlation within the [OI]5577 covariance 
group is that with the O2 band (r=0.85), followed by the blue 
bands Bl (r=0.61), N I 5200A (r=0.56), B2 (r=0.55), C3 
(r=0.43) and C2 (r=0.41). The correlation with the other con- 
tinuum regions is weaker (r <0.4). As for the Na I D group, be- 
sides the very tight correlations existing between the OH bands 
(r >0.97), the strongest correlation is observed between Na I D 
and OH(8-3) (r=0.75), followed by the other OH bands. Na I D 
correlates rather well with the continuum ranges C2 (r=0.58) 
and C6 (r=0.63) and with the O2 band (r=0.61). As anticipated 
in the previous section, the only meaningful correlation found 
in the third covariance group is between [OI]6300 and the N I 
5200 line (r=0.95). This is actually one of the tightest correla- 
tions found among all features. Finally, all continuum bands CI 
to C6 are well correlated with each other and with the blue bands 
B 1 to B4, which show as well a very strong mutual correlation. 

Even though the analysis is certainly hampered by the 
smaller sample, I have run a similar exploratory study for the 
blue emission features described in Sec. |6.3| The results are pre- 
sented in Table |9] The first interesting fact is that the blue emis- 
sion features bf 1 to bf4 are uncorrected with the CB continuum 
(|r| <0. 15), while some correlation is seen for bf5 (r ^0.4). On 
the other hand, all blue features show a significant mutual cor- 
relation, which ranges from a minimum (r ^0.4) for bf2 and 
bf5 to a maximum (r ^0.9) for bf3 and bf4. In general, bf5 is 
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Table 8. Linear correlation factors in the log F-log F plane. Only dark time data have been used. The last two columns and rows 
report the correlation with daily (SFD^/) and monthly (SFD,,,) solar flux density averages in the semi-log plane. 



CI 

0.15 
0.79 
0.14 
0.94 
0.31 
0.82 
0.02 
0.88 
0.23 
0.14 
0.15 
0.67 
0.14 
0.16 
0.87 
0.95 
0.99 
1.00 
0.60 



5200 

"OHT 

0.36 
0.43 
0.29 
0.56 
0.31 
0.95 
0.22 
0.34 
0.32 
0.32 
0.35 
0.29 
0.21 
0.31 
0.22 
0.16 
0.17 
0.60 



C2 
0.79 
0.36 

0.58 
0.94 
0.41 
0.99 
0.28 
0.96 
0.40 
0.32 
0.37 
0.92 
0.46 
0.35 
0.75 
0.77 
0.76 
0.79 
0.90 



NalD 
0.14 
0.43 
0.58 

0.35 
0.41 
0.52 
0.32 
0.41 
0.70 
0.72 
0.75 
0.63 
0.61 
0.61 
0.19 
0.14 
0.15 
0.16 
0.72 



C3 
0.94 
0.29 
0.94 
0.35 

0.43 
0.97 
0.17 
0.98 
0.21 
0.12 
0.16 
0.86 
0.33 
0.32 
0.85 
0.89 
0.92 
0.94 
0.80 



5577 
0.31 
0.56 
0.41 
0.41 
0.43 

0.38 
0.29 
0.31 
0.31 
0.29 
0.32 
0.26 
0.25 
0.85 
0.61 
0.55 
0.37 
0.36 
0.59 



C4 
0.82 
0.31 
0.99 
0.52 
0.97 
0.38 

0.26 
0.99 
0.36 
0.28 
0.32 
0.93 
0.44 
0.33 



6300 
0.02 
0.95 
0.28 
0.32 
0.17 
0.29 
0.26 

0.19 
0.15 
0.13 
0.13 
0.28 
0.15 
0.01 



C5 
0.88 
0.22 
0.96 
0.41 
0.98 
0.31 
0.99 
0.19 

0.27 
0.19 
0.23 
0.91 
0.37 
0.26 



(6-1) 
0.23 
0.34 
0.40 
0.70 
0.21 
0.31 
0.36 
0.15 
0.27 

0.98 
0.97 
0.67 
0.97 
0.53 



(7-2) 
0.14 
0.32 
0.32 
0.72 
0.12 
0.29 
0.28 
0.13 
0.19 
0.98 

1.00 
0.63 
0.99 
0.57 



(8-3) 
0.15 
0.32 
0.37 
0.75 
0.16 
0.32 
0.32 
0.13 
0.23 
0.97 
1.00 

0.65 
0.98 
0.60 



C6 
067" 
0.35 
0.92 
0.63 
0.86 
0.26 
0.93 
0.28 
0.91 
0.67 
0.63 
0.65 

0.64 
0.34 



(6-2) 
0.14 
0.29 
0.46 
0.61 
0.33 
0.25 
0.44 
0.15 
0.37 
0.97 
0.99 
0.98 
0.64 

0.47 



O2 
0.16 
0.21 
0.35 
0.61 
0.32 
0.85 
0.33 
0.01 
0.26 
0.53 
0.57 
0.60 
0.34 
0.47 



Bl 
OW 
0.31 
0.75 
0.19 
0.85 
0.61 



0.81 0.02 0.87 0.25 0.17 0.18 0.68 0.16 
0.83 0.02 0.88 0.25 0.17 0.18 0.68 0.16 
0.90 0.49 0.83 0.61 0.55 0.59 0.90 0.63 



- 0.96 
0.24 0.91 
0.22 0.89 
0.52 - 



B2 B3 B4 
0.95 0.99 1.00 
0.22 0.16 0.17 
0.77 0.76 0.79 
0.14 0.15 0.16 
0.89 0.92 0.94 
0.55 0.37 0.36 

- 0.81 0.83 

- 0.02 0.02 

- 0.87 0.88 

- 0.25 0.25 

- 0.17 0.17 

- 0.18 0.18 

- 0.68 0.68 

- 0.16 0.16 

- 0.24 0.22 
0.96 0.91 0.89 

0.99 0.97 
0.99 0.99 
0.97 0.99 

- 0.61 0.62 



Broad 
""060" 
0.60 
0.90 
0.72 
0.80 
0.59 
0.90 
0.49 
0.83 
0.61 
0.55 
0.59 
0.90 
0.63 
0.52 



0.61 
0.62 



SFDrf 
0.25 
0.31 
0.12 
0.03 
0.23 
0.30 
0.13 
0.09 
0.14 
0.02 
0.06 
0.04 
0.02 
0.06 
0.22 
0.28 
0.28 
0.24 
0.27 
0.14 



0.25 
0.23 



0.31 
0.45 



0.12 
0.14 



0.03 
0.07 



0.23 
0.22 



0.30 
0.49 



0.13 
0.13 



0.09 
0.23 



0.14 
0.10 



0.02 
0.07 



0.06 
0.04 



0.04 
0.06 



0.02 
0.05 



0.06 
0.10 



0.22 
0.41 



0.28 
0.25 



0.28 0.24 0.27 
0.24 0.25 0.26 



0.14 
0.24 



0.71 



o 
o 
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Fig. 23. Correlation between N 1 5200 and [01]5577 (upper 
panel) and [OI]6300 (lower panel). Only airmass correction has 
been applied. 



the feature that shows the weakest correlation with the remain- 
ing blue features. Additionally, it displays the strongest corre- 
lation with the [OI] 5577 feature (r ^0.8) and solar activity. 
Indeed, the correlation with [OI] 5577 is rather marked for all 
bf 's (r >0.5), suggesting that these features might belong to the 
green-line group of Barbier ( |1956| l. Nevertheless, bfl and bf2 
show a similarly marked correlation to the N 1 5200 line (r >0.5), 
that belongs to the red-line group. This suggests a partial corre- 



lation between the abrupt micro-auroral events undergone by the 
[OI]6300,6364A doublet and the blue features activity. 

As far as the correlation with solar activity is concerned, it is 
interesting to note that while the blue pseudo-continuum shows 
a very similar correlation factor with S FDj and S FD,„, for the 
blue features bf2 and bf3 this is significantly larger when the 
daily averages are used. On the contrary, bf5 shows a stronger 
correlation with the monthly averages. Because of the limited 
sample, though, these results have to be taken with some caution. 
In fact, for instance, the [OI]5577 line shows a suspiciously high 
correlation factor with S FDj (0.54), at variance with the value 
derived from the whole data set (0.30; see Table|8]l. 

8. Short timescale variations 

In general, the spectroscopic time coverage of the present sam- 
ple is sparse and this makes the study of short time scale vari- 
ations (minutes to hours) rather difficult. Interestingly, the sam- 
ple presented here includes some data sets for which the spec- 
troscopic observations went on for a significant fraction of the 
night. The most complete case is illustrated in Fig. 



24 which 



presents the behavior of some selected features during the night 
of 2001-04-22 as a function of time elapsed from the end of 
evening twilight, with an average interval of thirty minutes. This 
sequence definitely shows that the variations are smooth, with 
typical timescales of the order of several hours. As expected, the 
most marked change is seen in the [OI]6300 feature that, dur- 
ing the first hours of the night, went through one of the intensity 
enhancements typical for sites close to the geomagnetic equator. 

The smooth time evolution seen in the emission features 
matches the behavior observed in broad band photometry when 
long time series are available (see Paper I, Sec. 6). 

9. Discussion 

The analysis presented in this paper gives a clear picture of the 
complexity shown by the nightglow fluctuations, most of which 
remain unexplained. If it is well established that several features 
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CI 5200 C3 5577 CB bfl bf2 bf3 bf4 bf5 SFD^ SFD„, 



CI 
5200 

C3 

5577 

CB 

bfl 

bf2 

bO 

bf4 

bf5 



0.21 
0.97 
0.40 
1.00 
0.05 
0.01 
0.15 
0.15 
0.40 



0.21 

0.32 
0.59 
0.20 
0.56 
0.51 
0.21 
0.12 
0.13 



0.97 
0.32 

0.48 
0.96 
0.15 
0.19 
0.23 
0.07 
0.46 



0.40 
0.59 
0.48 

0.38 
0.65 
0.47 
0.54 
0.56 
0.80 



1.00 0.05 
0.20 0.56 
0.96 0.15 
0.38 0.65 
0.05 

0.05 

0.00 0.82 
0.15 0.76 
0.14 0.54 
0.38 0.50 



0.01 
0.51 
0.19 
0.47 
0.00 
0.82 

0.80 
0.69 
0.44 



0.15 
0.21 
0.23 
0.54 
0.15 
0.76 
0.80 

0.91 
0.62 



0.15 
0.12 
0.07 
0.56 
0.14 



0.40 
0.13 
0.46 
0.80 
0.38 



0.54 0.50 
0.69 0.44 
0.91 0.62 
0.69 

0.69 



0.32 
0.49 
0.38 
0.54 
0.30 
0.39 
0.32 
0.21 
0.26 
0.50 



SFDrf 
SFD,„ 



0.30 
0.46 
0.36 
0.45 
0.28 
0.30 
0.08 
0.02 
0.18 
0.32 



0.32 
0.30 



0.49 
0.46 



0.38 
0.36 



0.54 
0.45 



0.30 0.39 
0.28 0.30 



0.32 
0.08 



0.21 
0.02 



0.26 0.50 
0.18 0.32 



0.79 



0.79 



4.0 
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Fig. 24. Sequence of spectral measurements for some of the fea- 
tures studied in this work. The data were obtained with FORS 1 
on 2001-04-22 (G300V+GG435). Exposure times ranged from 
10 to 20 minutes. Several sky targets at different airmasses were 
observed. The vertical arrows indicate the times of target change. 
Airmasses spanned from 1.2 to 1.9. The vertical dashed line 
marks the start of morning astronomical twilight. The horizon- 
tal dashed lines mark the average value, while the dotted lines 
indicate +20% levels. 



show a definite correlation with solar activity, to which they react 
with timescales of the order of a couple of weeks, more thorough 
investigations need to be performed in order to better understand 
the link between space weather and the phenomena taking part 
in the upper layers of Earth's atmosphere. 

In all the studies of the nightglow in the astronomical con- 
text, the radio flux at 10.2 cm has been used as the only proxy 
for the solar activity. Nevertheless, other transient solar phenom- 
ena might have some impact on the night sky brightness, like 
Flares, Coronal Mass Ejections and Solar Proton Events (see 
Hanslmeier 12007 1 for a very recent review on solar phenomena). 
During these events, large amounts of energetic charged particles 
are released and, in the case they interact with the Earth's mag- 
netic field, they cause a series of geomagnetic effects, including 



aurorae. Therefore, charged particles are potential responsibles 
for at least some of the observed nightglow fluctuations, includ- 
ing the SAO. 

Even though this will require a dedicated analysis, I have 
run a preliminary study using the proton flux measured by the 
CELIAS/MTOF Proton Monitor (Ipavich et al. 1998^ on board 
of soHq5 orbiting on the Earth-Sun line at about LSxlO*" km 
away from Earth. The proton number flux (PNF) has been de- 
rived multiplying the proton density by the proton velocity mea- 
sured by CELIAS/MTOF and it is plotted in Fig. [25] for the rel- 
evant time interval (upper panel). 

The PNF average value during this time span is 2.6x10^ 
cm" s" (corresponding to a mass loss rate of -2x10"'"^ Mq 
yr"'), while the yearly averages range from 1.9 to 3.5x10** cm"^ 
s 7 The maximum average value was reached in July 2002, i.e. 
some time after the secondary maximum seen in the radio flux. 
The proton velocity ranges from 270 to 1000 km s"' (median 
value 440 km s"'), while the number density ranges from 0.1 
to 73 cm"-' (median value 4.7 cm"^). Looking in more detail at 
the SOHO data, one notices that a number of isolated and short 
duration peaks are present in the PNF. In fact, in some cases, 
values as high as 4x10^ cm"^ s"' are reached. Already a look at 
Fig. [25l (upper panel) shows that the frequency of these spikes is 
higher closer to the solar maximum. This can be seen more quan- 
titatively in the bottom panel of Fig. 25 where I have plotted the 
rate of what I will indicate as proton events (PEQ In this context, 
a PE has been defined as a group of adjacent CELIAS/MTOF 
measurements with PNF>7xlO^ cm"^ s"'. 

In principle, since the solar rotation axis is inclined by about 
7. 1 degrees with respect to the ecliptic (Balthasar, Stark & Wohl 
1987), and because the solar wind flux is known to be higher 
from the poles than the equator, a modulation of the solar wind 
as seen from Earth is indeed expected. More precisely, the pro- 
ton flux should be maximum around March 5 and September 
5 when Earth lies at its highest/lowest heliographic latitude, re- 
spectively. For the same reason, the proton flux should be min- 
imum around January 5 and June 5. This prediction can be eas- 
ily compared to the real data using the CELIAS/MTOF mea- 
surements and looking at their behavior as a function of time 
elapsed since the beginning of the year. The result is shown in 
Fig. 26 which was produced using data obtained towards solar 
minimum (January 2004 to January 2007), for a total of about 
28,000 data points. Indeed the PNF shows a SAO, with maxima 
in April and November, i.e. significantly shifted in time with re- 
spect to the epochs of maximum/minimum heliographic Earth's 



* Data can be downloaded from|http://umtof.umd.edu/pin/cm/| 
' Note that these events differ from the Solar Proton Events, which 
are bursts of relativistic protons, with energies larger than 10 MeV. 
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Fig. 25. Upper panel: CELIAS/MTOF proton number flux. 
Central panel: solar radio flux. Lower panel: CELIAS/MTOF 
proton events (PF>7xlO^ cnT^ s"' ). 



latitude. Remarkably, the SAO observed in the PNF appears to 
be in phase with the SAO detected both in the broad band data 
(Sec. [5]l and in the emission features (Sec. [6|l. Interestingly, a 
similar plot for the 10.2 cm radio flux does not show any clear 
trace of a SAO. Also, the evidence of a SAO in the PNF for the 
years around the maximum of solar cycle n. 23 becomes weaker 
A possible explanation is that the more frequent and probably 
more stochastic PEs tend to dominate over the smoother SAO 
when the sun is more active. 

Besides the smooth, long-term variation possibly produced 
by the modulation of the average proton fluence, it is reasonable 
to think that isolated PEs might be the cause of sporadic night- 
glow enhancements. To explore this possibility, 1 have run a cor- 
relation analysis similar to the one described in Sec. 4.1 for the 
solar radio flux. The results are in general rather noisy and show 
that the correlation with PNF is always low (|r| <0.3), both for 
short (Af=0.5 days) and long (Af=30 days) time windows. An 
example for the V passband is presented in Fig. [27](upper panel) 
for an averaging window Af=10 days. The correlation coefficient 
reaches a peak ~0.3 for t ^10 days, it decreases and it sud- 
denly drops for t >27 days. Even though the result is not really 
convincing, it might indicate that the response of the night sky 
brightness to PEs takes place with some time delay and that what 
matters is the proton flux behavior during the last solar rotation. 
It must be noticed that the PNF shows a rather marked recur- 
rence, with a period of 9.1 days, that coincides with one third of 
the solar rotation period (see Fig. 27 lower panel). Similarly to 
what happens for the radio data (see Sec. |4. this probably cre- 
ates spurious correlation peaks. Interestingly, shorter-timescale 
periodic variations in the geomagnetic activity have been de- 
tected for 13.3, 9.1 and 6.9 days (see for instance Hauska, Abdel- 
Wahab & Dyring 197311. 

Since the typical speed of solar wind is ~450 km s"', the 
swarms of particles released during the PEs reach the Earth 
about one hour after being detected by SOHO. Therefore, the 



Fig. 26. Seasonal variation of the MTOF/CELIAS proton num- 
ber flux. Only data from 2004-01-01 to 2007-01-01 are plotted. 
The filled circles trace the monthly averages. No correction for 
the 1 1 -years solar cycle has been apphed. 

delay in the reaction is completely due to processes taking place 
within the Earth's atmosphere. 

As a last check for intermediate timescale periodic varia- 
tions, I have investigated the correlation with the Moon mo- 
tion. In fact, it has been suggested that atmospheric tides might 
induce recurrent oscillations in the night sky brightness (see 
Chamberlain 1961 and references therein). For this purpose, I 
have investigated possible links between the broad band, dark 
time measurements and lunar phase or lunar hour angle. No sig- 
nificant correlation has been found, in agreement with a similar 
analysis run by Mattila et al. d 19961 1. 

10. Conclusions 

In this paper I have presented a photometric and spectroscopic 
analysis of the optical night sky emission at Cerro Paranal in the 
time interval April 2001 - January 2007. The main conclusions 
of this work can be summarized as follows; 

- The UB VRI night sky brightness is well correlated with solar 
activity. The correlation is maximum in U and minimum in 
/. 

- The excursion between sunspot maximum and minimum of 
solar cycle n. 23 is 0.6, 0.3, 0.3, 0.4 and 0.2 mag arcsec"^ in 
U, B, V, R and / passbands, respectively. 

- There are indications that the effects of solar activity on 
nightglow are not identical across different solar cycles. 

- The reaction of the nightglow to the solar variations seems 
to take place with a time delay of the order of 2-3 weeks. 

- The R night sky brightness seems to react with a much 
shorter delay, of the order of a few days. 

- V,R and / measurements show a clear semi-annual oscilla- 
tion (SAO), with a typical peak-to-peak variation of ~0.5 
mag arcsec"^. For the B passband this oscillation is, if 
present, much smaller. 
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Fig. 27. Upper panel: Linear correlation coefficient between dark 
time V sky brightness and PNF as a function of time delay r, 
computed for At=10 days. Lower panel: example of PNF peri- 
odicity. The dotted curve is a replica of the original data shifted 
by 9.1 days. 

- Maxima and minima of the SAO are out of phase with re- 
spect to the Equinoxes and Solstices. 

- All main emission features show a SAO, very similar to the 
well known seasonal oscillation of the Na 1 D doublet. 

- [01] 5577 and NI 5200A show the strongest correlation with 
solar activity. For [OI] 5577, the maximum correlation is 
found for a time delay of 15 days. 

- [01] 6300Aand NI 5200A show a very tight mutual corre- 
lation. Nevertheless, [OI]6300 displays a weaker correlation 
with solar activity. 

- Flux variations in the OH bands are very strongly correlated 
with each other and do not show any correlation with solar 
activity. 

- The main emission features, both atomic and molecular, 
show smooth flux variations on time scales of hours. 

- A preliminary and exploratory analysis of the possible con- 
nection between nightglow and flux of charged particles 
from the Sun has shown that there is a weak correlation. 

- The night sky emission seems to react, with a delay of about 
10 days, to variations in the proton number flux. 

- No correlation is found between the dark time, broad band 
night sky brightness and moon phase or moon hour angle. 

Future investigations, with even larger databases, will have 
to address the possible relations with other solar phenomena, like 
Coronal Mass Ejections, Flares and Solar Proton Events, in an 
attempt to connect the observed short timescale variations of the 
nightglow with space weather. 
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